16.451 Lecture 17: Beta Decay Spectrum {Ref: F&H ch 11 ) Nov. 3, 2005 


n — » p + e + v e 

(and related processes...) 


Goals: 


• understand the shape of the energy spectrum 

• total decay rate sheds light on the underlying 

weak interaction mechanism 



Starting point: "Fermi's Golden Rule" again! (lecture 6) 

matrix element 
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M if =J Yf V{r)Yi d r 

4 
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Pf = dn/dEf 


(transitions / sec) 


density of states 


Matrix element M lf ? 


Simplest model is to take a pointlike interaction with an overall energy scale "G": 

(Fermi, 1934 - almost right!) 



M if = G | f e {r) <//„,,■ (r) d 3 r 


(the interaction is proportional to the wavefunction 
overlap of initial and final state particles at the same 
point in space) 


P 

e 

w- y e 

Standard Model description: an extended 
interaction, but the range is only about 0.002 fm 
which is just about zero! 

n 


The Standard Model can ’predict' the value of G 
in terms of model parameters, whereas in Fermi's 
theory, it remains to be determined from experiment. 


i 






Spin considerations: electron and neutrino 3 

There are two possibilities for the angular momentum coupling of the two leptons! 

s e +s v = S tot) => S = 0 or 1 

For neutron decay: n —> p + e + V e 

Angular momentum: i = i + S both can contribute to neutron decay! 

Subtle point: because the leptons are emitted with a definite helicity 

we can deduce a correlation between their directions of motion in the two cases: 


Fermi Decay, S = 0: 


Samow- Teller decay, S = 1: 

s e ( LH ) Pe 


S e (LH) Pe 




S v ( RH ) ^ 


Sy (RH) j: 

Hv 

e- and v travel in the same direction 


e- and v travel in opposite directions 


More on spin correlations: 4 

Fermi Decay, S = 0: n — » p + e~ + V e 

Leptons travel in the same direction:: 

p P 

Recoiling proton spin is in the same direction 
as the initial neutron spin. 

S = 0 

Gamow - Teller Decay, S = 1: 

Leptons travel in the opposite direction: 


S = 1 



Recoiling proton spin is in the opposite direction 
as the initial neutron spin, i.e. a “spin flip" 
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How to proceed? 5 


As before, assume a pointlike interaction, but allow for different coupling constants 
for the Fermi (F) and Gamow-Teller (GT) cases. 


Fermi case, S = 0: (coupling constant: "G v " because the operator transforms 

like a space vector.) 

M if = Mp = Gy jVp,f(r )fe{r)<t>l (r )y/ ni (r) d 3 r 

Gamow-Teller, S = 1: (coupling constant: "G A " because the operator transforms 

like an axial vector, i.e. like angular momentum.) 

M if = m gt = g a \v*p, f (j) f e {r) <j)* v {r)\i/ ni {r) d 3 r 

Experimentally, the coupling constants are very similar: 

Ga / Gy = -1.25 

(These are evaluated by comparing different nuclear beta decay transitions where 
angular momentum conservation restricts the total lepton spin states that contribute) 


Transition rate: 


for the neutron: 


n ^ p + e + v e (can proceed with S = 0 or 1) 

A if = ^ \M if \ 2 Pf ~ (Gy +3 G 2 a ) 
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since there are 3 times as many ways for the leptons to be emitted with S = 1 
(m s = 1, 0, -1) as with S = 0. 

For now, let us work out a generic matrix element, since the expressions are the 
same for both apart from the coupling constants: 

M if = G jV* j(r)f e (r)f v (r)y/ ni (r) d 3 r 



Q neutrino 


Electron and neutrino are represented by 
plane wave functions of definite momentum: 


<t>e(f) 


e ip*r/h 


etc. 
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Matrix element: 


7 


Pr 

o 

recoil 



M if = G J i/ pJ (r) </>*,(?) f v (r) y ni (r) d 3 r 


Pr = ~(P + Q) 


Q neutrino 



l) e -i(p+q)»f/h 



e i P R *r/h 


The integral for M lf extends over all space regions for which the nucleon wave functions 
(n,p) are non-zero: R max ~ 1 fm (in nuclei, use R - 1.2 A 1/3 fm ) ... 

But, the recoil momentum p R is no larger than the Q-value for the reaction, - MeV ... 


p R *r/fi < 


QK 


he 


1 MeV.fm 

« 1 

197 MeV.fm 


l=> 



This is a great simplification: the lepton wave functions are just a constant over the 
region of space that matters to calculate the matrix element © 


Summary so far for the transition rate calculation: 
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The remaining integral is referred to as the nuclear matrix element: 


M 


nuclear 


= jVf,p(r) y/ Un (r)d 3 r 


When beta decay occurs in a nucleus, the initial and final wave functions of the proton 
and neutron need not be exactly the same, so in general: 


M 


nuclear 


< 1 


However, in the case of the free neutron, there are no complicated nuclear structure 
effects, and so the matrix element is identically 1: 


When this occurs in a nucleus, the beta decay 
rate is the fastest possible, and the 
transition is classified as "superallowed" 


A if - 


2 n_ G_ 

h y 2 


Pf 
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Density of states factor: 


Just like the calculation we did for electron scattering, but now there are two light 
particles in the final state! 


Pr 


recoil 


o 


electron 


We want to work out 
the number of equivalent 
final states within energy 
interval dE f of E f . 


pf 


dn 

dE f 


4 neutrino 

Final state momenta are quantized in volume V (lecture 6) 


dn = dn e x dn v = I An p 2 dp x \ An q 2 dq V 


E j — Me + me + K R + K e + cq 

But the nucleon is much heavier than the other particles: K R = (p R ) 2 / 2 M = 0 


A dEf\ = c dq 

~ r 1 I K p =const 


Put this all together for the density of states factor: 
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Pf - 


dn 

dE f 


dn e x dn v 


c dq 


= (4^) 2 ^ dp 
he 


And finally, for the transition rate: 



Exercise: plug in all the units and check that the transition rate is in sec 1 


Notice: This is actually a partial decay rate, because the electron momentum p is 
specif ied explicitly. X if here gives the rate at which the decay occurs for a given 
electron momentum falling within dp of p -» this predicts the momentum spectrum! 
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Electron momentum spectrum: 


N(p)dp = N 0 Ajf = (const.) x p 2 q 2 dp 
=^> N(p) = (const.) x p 2 q 2 


approx: K R = 0 

. / 

= (const) xp 2 ( Q-K e ) 2 



' o % 

recoil 


Q neutrino 


Predicted spectral shapes, 
Krane, figure 9.2: 

(plotted for Q = 2.5 MeV, 
not the neutron!) 

(Note: max. K(e') = Q ) 



Comparison to Reality: (Krane, fig. 9.3: e+ and e- decays of 64 Cu) 12 


Too many low energy e 



Too few low energy e + 
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Discrepancy: neglect of Coulomb effects in the final state. 


Key point: Coulomb distortions of the energy spectra occur AFTER the electron/positron 
are emitted in the weak decay process. 

Modif ied density of electron/positron states: 



"Fermi function", depends on 
the charge Z' of the "daughter 
nucleus" (final state) and the 
electron/positron momentum 


Approximate correction factor for p± decay: 




x 

1- e~ x ’ 


x 


+ 


2 n a Z' 

P 


p 


v e 2 1 

— , a- 

c 4ns Pic 137 


Modified electron/positron spectrum prediction: 


N(p) = Cp 2 ( Q-K e f F ± (Z', p) , C = 


: 3 


M 


nucl 


Fermi-Kurie Plot 


14 


Idea: for "allowed decays", corresponding to our approximation: e iPR-rlh = x 

inside the nucleus, the electron energy spectrum can be "linearized" if one 
accounts for the Coulomb distortion via the Fermi f unction F(Z',p): 



66 Ga -> 66 Zn decoy \ Phys. Rev. 129, 1782 


N(p) 

V p 2 f(z\p) 


Q -K e 


linear function, 
endpoint Q 


K e + m e 
m e 
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Neutrino Mass effect: 3 H — » 3 He + e + v e (Q = 18.6 keV) 


Idea: shape of the electron energy spectrum near the endpoint (Q) is sensitive to the 
mass of the electron antineutrino: 


Pr 

recoil 


electron 



neutrino 


recall: Q = K R + K e + K y 

When K e = Q, K R = K v ^ 0. if m v * 0, then 
in this limit, mass effects are most pronounced. 



energy E [keV] 


E-E g [eV] 


Figure 2: The electron energy spectrum of tritium ft decay: (a) complete and (b) narrow region 
around endpoint Eq. The ft spectrum is shown for neutrino masses of 0 and 1 eV. 


f igure from "KATRIN" proposal, 2001: (Karlsruhe Tritium Neutrino expt.) 
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Figure 4: Results of tritium ft decay experiments on the observable ni^ over the last decade. 

• Best direct upper limit: m v < 2.2 eV 

• from Sudbury neutrino observatory and other experiments, we have convincing 

indirect evidence of nonzero neutrino mass that is much smaller than this 



Web sites: http://cupp.oulu.f i/neutrino/nd-mass.html , http://www-ik.fzk.de/~katrin/ 
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